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Abstract 

We review some known facts in the transport theory of mesoscopic systems, including 
counting statistics, and discuss its relation with the mathematical treatment of open 
systems. 

1 Introduction 

The aim of these notes is to introduce to some theoretical developments concerning transport 
in mesoscopic systems. More specifically, we intend to show how concepts and tools from 
mathematical physics provide ways and means to put some recent, fundamental results on 
counting statistics on rigorous ground and in a natural setting. We will draw on concepts 
like C*-algebras, which have been often used in the mathematical treatment of systems out of 
equilibrium, see e.g. [7], but also on tools like Fredholm determinants, which have been used 
for renormalization purposes in quantum field theory. Before going into mathematical details 
we will review some of the more familiar aspects of transport, and notably noise. That will 
provide some examples on which to later illustrate the theory. 

These notes are not intended for the expert. On the contrary, the style might be overly 
pedagogical. 

2 Noises 

Consider two leads joined by a resistor. The value of its conductance, G, is to be meant, 
for the sake of precision, as corresponding to a two-terminal arrangement, meaning that the 
voltage V is identified with the difference of chemical potentials between right movers on 
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the left of the resistor and left movers on its right. We are interested in the average charge 
{Q) transported across the resistance in a time T, and in the variance ((Q^)) = (Q^) — {Q)"^-, 
equivalently in the current (Q) /T and in the noise ((Q^))/T. 
There are two types of noises: 

i. Equilibrium, or thermal, noise occurs in the absence of voltage, V = 0, and at positive 
temperature > 0. Then 

{Q) = o, ^ = 1^- (1) 

(Johnson [8], Nyquist [H]). This is an early instance of the fluctuation-dissipation 
theorem, those words being here represented as noise and conductance. 

ii. Non-equilibrium, or shot, noise occurs in the reverse situation: V ^ 0, I3~^ = 0. Ohm's 
law states {Q)/T = GV, while for the noise different expressions (corresponding to 
different situations) are available: (a) classical shot noise 

m = e{Q) (2) 

(Schottky [21]), where e is the charge of the carriers, say electrons. The result is inter- 
preted on the basis of the Poisson distribution 

A" 

Pn = e-^—, (n = 0,1,2,...) 
n\ 

of parameter A, for which 

(n) = A, {{n'))=\. 

Assuming that electrons arrive independently of one another, the number n of electrons 
collected in time T is so distributed, whence ([2]) for Q = en. 

(b) quantum shot noise. Consider the leads and the resistor as modelled by a 1- 
dimensional scattering problem with matrix 

s=(;f,). (3) 

where r, t (resp. r', t') are the reflection and transmission amplitudes from the left and 
from the right. Then 

)=e(g)(l-|tp) (4) 
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(Khlus [in], Lesovik [12] )• In this case the result may be attributed to a binomial 
distribution with the success probability p and with attempts: 



This yields (jl]) for p = |tp being the probability of transmission. For small p it reduces 
to ([2]). It should be noticed that in the case of thermal noise, the origin of fluctuations 
is in the source of electrons, or in the incoming flow, depending on the point of view. By 
contrast, in the interpretation of the quantum shot noise the flow is assumed ordered, 
as signified by the fixed number of attempts, and fluctuations arise only because of the 
uncertainty of transmission. 

We refer to [6] for a more complete exposition of these matters. We conclude the section 
by recalling that noises are quantitative evidence to atomism. Thermal noise determines 
l3~^ = fc-temperature, and hence Boltzmann's constant k as well as Avogadro's number Nq = 
R/k (somewhat in analogy to its determination from Brownian motion [23l [5]). Shot noise 
determines the charge of carriers. In some instances of the fractional quantum Hall effect this 
yielded e/3 [I9] or e/5 [H]. 

3 A setup for counting statistics 

Before engaging in quantum mechanical computations of the transported charge we should 
describe how it is measured, at least in the sense of a thought experiment. Consider a device 
(dot, resistor, or the like) connected to several leads, or reservoirs, one of which is distinguished 
('the lead'). The measurement protocol consists of three steps: 

• measure the charge present initially in the lead, given a prepared state of the whole 



• act on the system during some time by driving its controls (like gate voltages in a dot), 
but not by performing measurements. This includes the possibility of just waiting. 

• measure the charge present in the lead finally. 

The transported charge is then identified as the difference, n, of the outcome of the measure- 
ments. For simplicity we assume that n takes only integer values, interpreted as the number 




(n) = Np , {{n^)) = Np{l - p) . 



system. 
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of transferred electrons. Let p„ be the corresponding probabilities. They are conveniently 
encoded in the generating function 



of the moments of the distribution, 

in") 



. d 



A=0 



Similarly, logx(A) generates the cumulants {{n'')), inductively defined by 



in 



(5) 



where V = {ai, . . . , a^} runs over all partitions of {1, . . . , k}. Alternative protocols with 
measurements extending over time will be discussed later. 



4 Quantum description 

The three steps of the procedure just described can easily be implemented quantum mechan- 
ically by means of two projective measurements and by a Hamiltonian evolution in between. 

Let H be the Hilbert space of pure states of a system, p a density matrix representing 
a mixed state, and A = ^2 - aiPi an observable with its spectral decomposition. A single 
measurement of A is associated, at least practically, with the 'collapse of the wave function' 
resulting in the replacement 

i 

where ii{PipPi) = ii{pPi) is the probability for the outcome a^. Two measurements of A, 
separated by an evolution given as a unitary U , result in the replacement [22] 

p - J]P,f/P.pP.f/*P,-, (7) 
hi 

where ti {PjU Pi pPiU*Pj) = tiiU* PjU PipPi) is the probability of the history (aj,aj) of out- 
comes. We can so compute the moment generating function ([5]): 

;^(A) = ^tr(t/*P,t/PipPi)e'^^"'""'^ = Y^iAU*e'^^UP^pPi)e-'^''^ . (8) 

i,3 i 
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The expression simplifies if 



[Ap] =0; 



(9) 



tfien PipPi = PiP, wlience tlie r.li.s. of ([6]) still equals p (no collapse at first measurement) 
and 



5 Independent, uncorrelated fermions 

We intend to apply ffTOj) to many-body systems consisting of fermionic particles which are 
uncorrelated in the initial state. The particles shall contribute additively to the observable 
to be considered and evolve independently of one another. The ingredients can therefore be 
specified at the level of a single particle. At the risk of confusion we denote them like the 
related objects in the previous section: A Hilbert space H with operators A, U, p. However, the 
meaning of p is now that of a 1-particle density matrix < p < 1 specifying an uncorrelated 
many-particle state, to the extent permitted by the Pauli principle: any eigenstate of \i>) of 
p, p\iy) = I'll'), is occupied in the many-particle state with probability given by its eigenvalue 
u. Common examples are the vacuum p = and, in terms of a single-particle Hamiltonian 
H, the Fermi-Dirac distribution p = (1 + e^^)~^ or its zero temperature limit, 0, the 

Fermi sea p = Q{—H). 

The corresponding many-particle objects are obtained through second quantization, which 
amounts to the following replacements: 




tr(t/*e^^^f/e-^^^p) . 



(10) 



X(A) 



(11) 



oo 



n 




(Fock space) 



(12) 



n=0 



A ^ dT{A) 
U T{U) 



(13) 
(14) 



where dr(yl) and T{U) act on the subspaces /\^ H C J-'(H) as 



n 




i=l 



T{U) =U^---®U. 
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Moreover, the state is replaced as 

P ^ , (P'=l-P). (15) 

Tr^(H)r(p/p') 

Indeed, if p splits with respect to = Tii © 7^2, then the many-body state f|T5|) factorizes 
w.r.t. jF(7i) = jF(7ii) ® ^(^2). In particular if p|i^) = z/|i^), this entails the following state 

on^[|z.)] = ©i=oA"[k)] 

. , _^ = z/lo + i^li, 

confirming that u is the occupation number of |z/). We note that 

Tr^(H)r(M) = det„(l + M), 

provided that M is a trace-class operator on 7i, in which case the r.h.s. is a Fredholm 
determinant. We will comment on that condition later. Under the replacements (11211151) 
the assumption [A, p] = is inherited by the corresponding second quantized observables, 
[dr{A),r{p/ p')] = 0. As a result fllOl) applies and becomes the Levitov-Lesovik formula 

X(A) = det(p' + e'^^*^^e-^^^p) . (16) 

Indeed, 

_ Tr^(H)(r(?7)-e'^'^^(-^)r(^)e-^^'^r(A)r(^/^.)) 

^ Tr^wr(p/p') 
_ TT^^H}^{U*e'^^Ue~'^'^p/p') _ detjl + U*e'^^Ue-'^^{p/ p')) 
~ Tr^(„)r(p/p') " det(l + (p/p')) 

= det(p' + f/*e'^^t/e-'^^p). 

Before discussing the mathematical fine points of (fT6|) . let us compute the first two cumu- 
lants of charge transport. In line with the discussion in the previous section, let A = Q be 
the projection onto single-particle states located in the distinguished lead. Then (ITBi) yields 

(Q) = -ix'(O) = trp(AQ) , 
HQ')) = -(logx)"(0) = trp(Ag)(l - p)AQ 

= tr(p(l - p){AQf) + ^tr(i[Ag, p])^ , (17) 



where AQ = U*QU — Q is the operator of transmitted charge. The split ( 1171) of the noise 
{{Q')) into two separately non negative contributions is of some interest ([3] by a different 
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approach, [T]): The commutator [A(5,p] expresses the uncertainty of transmission AQ in the 
given state p; the second term in (fTTI) may thus be viewed as shot noise. The factor p(l — p) 
expresses the fluctuation v{l — u) in the occupation of single particle states It refers to 
the initial state, or source, and its term may be viewed as thermal noise; indeed it vanishes 
for pure states, p = p^, while for p = (1 + e^^)~^ the energy width of p(l — p) is proportional 

to p-\ cf. m- 



6 Alternative approaches 

We present alternatives and variants of the two-step measurement procedures discussed in 
Sect. [31 We discuss them in the first quantized setting of Sect. HI The corresponding second 
quantized versions can then easily obtained from the replacements ( IT2lfT5|) . 

i) [11] One could envisage a single measurement of the difference U*AU — A. On the basis 
of ([7]) its generating function is 

X(A) = tr(e'^(^*^^-^)p) . 

It remains unclear how to realize a von Neumann measurement for this observable, since its 
two pieces are associated with two different times. Moreover, its second quantized version 

X(A) = det(p' + e'^(^*^^-^V) 
generates cumulants which, as a rule beginning with n = 3, differ from those of (fT6l) . 

ii) [2U] We keep the two- measurement setup, but refrain from making assumption i.e. 
, the first measurement is allowed to induce a "collapse of the wave function" . We do however 
assume that the eigenvalues of A are integers, in line with the application made at the end 
of the previous section, where A dr{Q) with Q a projection. Then ([8]) yields 

X(A) = Y,^l{U*e'^^UPnpPm)Smne-'''' 



= — I drtr(f/V^^f/e-^(^+")V'"^) 
2vr Jo 

by using 5mn = (27r)-i J^" dr e'^^'"""). 

iii) [13] Here neither ([5]) nor ([9]) is assumed. The system is coupled to a spin-| resulting 
in a total state space Ti. ® C^. Specifically, the total Hamiltonian is obtained by conjugating 
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the system Hamiltonian by e i2^'^°'3^ where A is a couphng constant and as a Pauh matrix; 
equivalently, the same holds true for the evolution U, which becomes 

[/ = e-'^^®"3(?7® l)e^^^®"^ 

We note that 

U{ij^\a)) = {U^.x^)®\a), (a = ±l), 

where as\a) = a\a) and Ux = e~^^^Ue^^^. The joint initial state is assumed of the form p®Pi 
with p being that of the system and 

that of the spin. The final state is U{p ® Pi)U* and, after tracing out the system, 

pi = tTnU{p(S) Pi)U* 

with matrix elements 

{a\pi\a')=tTiU^xpU:,x){(T\pM')- 

In other words, 

^ _ f P++ p+~xW\ 
V-+X(-A) P- ) 

with 

X(A) = ti{S^U*e-''^^uS'^p) . 

We remark that x(A) agrees with ffTOl) under the assumption ([9]) of the latter. It can be 
determined from the average spin precession, since {a\pi\a') reflects that measurement. On 
the other hand no probability interpretation, c/. ([5]), is available for x(A), since its Fourier 
transform is non-positive in general [9]. 

7 The thermodynamic hmit 

The derivation of fll6p was heuristic. It therefore seems appropriate to investigate whether 
the resulting determinant, cast as det(l + M), is well-defined, which is the case if M is a trace- 
class operator. This happens to be the case if the leads are of finite extent and the energy 
range finite, essentially because the single-particle Hilbert space becomes finite dimensional. 
While these conditions may be regarded as effectively met in practice, it is nevertheless useful 
to idealize these quantities as being infinite. There are two physical reasons for that. First, 
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any bound on these quantities ought to be irrelevant, because the transport occurs across the 
dot (compact in space) and near the Fermi energy (compact in energy); second, the infinite 
settings allows to conveniently formulate non-equilibrium stationary states. However this 
idealization needs some care. In fact, in the attempt of extending eq. f|T6|) to infinite systems, 
the determinant becomes ambiguous and ill-defined. The cure is a regularization which rests 
on the heuristic identity 

tT{U*pQU - pQ) = 0, (18) 

obtained by splitting the trace and using its cyclicity. It consists in multiplying the determi- 
nant by 

det{e-'^^'P^^) ■ det(e'^^^) = Q-^^^^iu-pQU-pQ) ^ i ^ (19) 

thereby placing one factor on each of its sides. The straightforward result is ([2], and in the 
zero-temperature case [16]) 

X(A) = detie-'^P^^'^'p'e'^P^ + e^HjQu p^-iXp'Q^ ^ ^20) 
where p' = 1 — p, pu = U*pU, and similarly for p^ and Qu- 

Remarks. 1. We observe a manifest particle-hole symmetry: 

Xp(A) = Xp'(-A) . 

2. We will see that the determinant (!20!) is Fredholm under reasonable hypotheses. 

3. The regularization bears some resemblance to det2(l + M) = det(l + M)e~^'^^^ , though the 
latter typically changes the value of the determinant. 

To the extent that the regularization is regarded as a modification at all, it affects only 
the first cumulant, because the term —iXti^puQu — pQ), which by f|T9l) has been added to the 
generating function logx(A), is linear in A. The mean is thus changed from {n) = ii p{Qu — Q) 
to (n) = tr(p — pu)Qu- In line with Sections |3] and O we interpret Q as the projection onto 
single-particle states in the distinguished lead and U as the evolution preserving the initial 
state p, except for changes in the dot. We then expect that Qu — <5 is non-trivial on states 
of any energy, while p — pu is so only on states which are located near the dot and near the 
Fermi energy. As a result, the second expression for (n), but not the first one, appears to be 
well-defined. 

8 A more basic approach 

The regularization ffTSjl remains an ad hoc procedure, though it may be motivated as a 
cancellation between right and left movers, see [2]. The point we wish to make here is that 
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eq. fl20p is obtained without any recourse to regularization if the second quantization is based 
upon a state of positive density (rather than the vacuum, cf. Sect. E]), as it is appropriate for 
an open system. 

To this end let us briefly recall the defining elements of quantum mechanics of infinitely 
many degrees of freedom: (local) observables are represented by elements of a C*-algebra A 
and states by normalized, positive, continuous linear functionals on A. A state u, together 
with its local perturbations, may be given a Hilbert space realization through the GNS con- 
struction: it consists of a Hilbert space H^j, a representation vr^ of A on H^j, and a cyclic 
vector G H^j such that 

u{A) = {Q^,7T^{A)Q^), [AeA). 

Notice that the state u is realized as a vector, Q^^, regardless of whether it is pure. Rather, 
it is pure iff the commutant ni^{A)' C CiHuj) is trivial. The closure of ni^{A) yields the von 
Neumann algebra 7i^{A). Besides of local observables n^{A) it also contains some global ones, 
whose existence and meaning presupposes uj. An example occurring in the following is the 
charge present in the (infinite) lead in excess of the (infinite) charge attributed to u. 

The C*-algebra of the problem at hand is A{Ti.), the algebra of canonical anti-commutation 
relations over the single-particle Hilbert space Ti. It is the algebra with unity generated by 
the elements a(/), a{f)* (anti-linear, resp. linear in / G 7i) satisfying 

{a(/), a*{g)} = (/, ^7)! , {a(/), a{g)} = = {a*(/), a*(^?)} . 

A unitary U induces a *-automorphism of the algebra by a{f) ^ ci{U f) (Bogoliubov auto- 
morphism). A single-particle density matrix < p < 1 defines a state uj on A{T-L) through 

uj{a*{f)a{g)) = {g, pf) , uj{a{f)a{g)) = = uia*if)a*ig)) 

and the use of Wick's lemma for the ccomputation of higher order correlators. States of this 
form are known as gauge-invariant quasi-free states; they describe uncorrelated fermions. It 
is possible to give an explicit construction of their GNS representation, known as Araki-Wyss 
representation, but we will not need it. 

For clarity we formulate the main result first for pure state and then for mixed states. In 
both cases we assume [p,Q] = 0, cf. IQ. 

Theorem (Pure states). Let p = p^. We assume that 

p-UpU* (21) 

is trace class. Then 
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i. The algebra automorphisms a(/) ^ ci{Uf) and a(/) a(e''^'^/) are unitarily imple- 
mentable: There exists (non-unique) unitaries U and e'^*^ on such that 



ii. Q is an observable, in the sense that any bounded function thereof is in 7ri^{A{Ti.)) . 

iii. The above properties define U uniquely up to a phase and Q up to an additive constant. 

iv. The generating function of cumulants, cf. ( ITTi) . equals the regularized determinant ( 120|) : 

{Q^, U*e'^^Ue-''^^Q^) = det(e-^^^^'3^p'e^^^^ + e^Vc/Qc/pe-iVQ) ^ 
where the determinant is Fredholm. 

Eq. (pri) demands that the evolution U preserves p, except for creating excitations of 
finite energy within an essentially finite region of space. This assumption is appropriate for 
the evolution induced by a compact device operating smoothly during a finite time interval. 

The generalization to mixed states is as follows. 

Theorem (Mixed states). Let < p < 1. Assume, instead of ([21]), that p^/^ - Up^^'^U* and 
(p')i/2 _ Ul^p'y/'^lJ* are trace class; moreover that 

ipp'Y^'Q (22) 

is, too. Then the above results (i-iv) hold true, upon replacing (iii) by 

• Properties (i-ii) define U uniquely up to left multiplication with an element from the 
commutant tTuj{-^{T^))\ and Q up to an additive constant. In particular, U*e^^^Ue~^^'^ 
is unaffected by the ambiguities. 

Notice that the most general case, < p < 1, is not covered. The physical origin of the 
extra assumption (l22l) needed in the mixed state case is as follows. In both cases, pure or 
mixed, the expected charge contained in a portion of the lead is of order of its length L, or zero 
if renormalized by subtraction of a background charge. In the pure case however, the Fermi sea 
is an eigenvector of the charge operator, while for the mixed state, the variance of the charge 
must itself be of order L, because the occupation of the one-particle states is fluctuating, due 
to pp' 7^ 0. Hence, in this latter situation, the measurement of the renormalized charge yields 
finite values only as long as L is finite, of which eq. fl22|) is a mathematical abstraction. In the 
limit L — > oo all but a finite part of the fluctuation of the source is affecting the transmitted 
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noise. That suggests perhaps that there is a better formulation of the result. Indeed, the 
expression for the transmitted noise, cf. the first term (IT7I) . is finite if {U*QU — Q){pp'y^'^ 
is trace class. This condition turns out to be sufficient for property (i), for making the 
determinant Fredholm and U*QU — Q an observable, but not for (ii, iii). 
For proofs we refer to [2]. 

9 An application 

We discuss a very simple application to illustrate the working of the regularization. The 
system consists of two leads in guise of circles of length T, joined at one point. Particles run 
in the positive sense along the circles C at velocity 1, whence it takes them time T to make a 
turn, and may scatter from one to the other circle at the junction. Initially states in the two 
circles are populated up to Fermi energies fiL < I^r- This is formalized as follows. The single 
particle Hilbert space is 



1-1 = L\C)®L\C) 3i) = 




the evolution over time T is 

{UiIj){x) = Sipix) 

with S as in ([3]). The momentum operator is p = —id/dx and the initial state is p = pi® Pr 
with Pi = 9{fii — p), {i = L,R). The projection onto the right lead is Q = © 1. 
Quite generally, for p = p^ a. pure state, eq. fl2Ul) reads 

X(A) = det«(l + (e~^" - l)QuPup' + (e^' - l)Qup'uP) , 
and in the present situation that determinant reduces to 

x{\)=detL2ic){l + {^'^-l)p'RPL\t?). 

It is to be noted that p'^Ph selects a finite energy interval, {fiLiP'R], unlike the determinant 
without regularization. Using eigenstates of momentum p G (27r/T)Z we find 

X(A)= n (l + (e^'-l)I^P) 

= {l-\t\^ + e'^\t\Y = {q + e%f 

with = ^{p \ Pl < P Pr} — ifJ'L — fJ'R)T/27r. This is a binomial distribution with 
probability p = |tp and attempts, reproducing (j4]). In particular, it yields Ohm's law 
(Q)/T = GV with G = (27r)-i|t|2 [THIl]. 
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